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We study the non-equilibrium version of the fluctuation-dissipation theorem (FDT) within the glass 
phase of Bouchaud's trap model. We incorporate into the model an arbitrary observable m and 
obtain its correlation and response functions in closed form. A limiting non-equilibrium FDT plot (of 
correlator vs response) is approached at long times for most choices of m, with energy-temperature 
FDT a notable exception. In contrast to standard mean field models, however, the shape of the 
plot depends nontrivially on the observable, and its slope varies continuously even though there 
is a single scaling of relaxation times with age. Non-equilibrium FDT plots can therefore not be 
used to define a meaningful effective temperature T^s in this model. Consequences for the wider 
applicability of an FDT-derived Toff are discussed. 



PACS: 05.20.-y; 05.40.-a; 05.70.Ln; 64.70.Pf 

Glassy systems relax extremely slowly at low temper- 
atures. They therefore remain far from equilibrium on 
very long time scales, and exhibit ageing [y: the time 
scale of response to an external perturbation increases 
with the age (time since preparation) i^ of the system. 
Time translational invariance (TTI) and the equilibrium 
fluctuation-dissipation theorem |g] (FDT) relating corre- 
lation and response functions break down. 

Let C(i,iw) = (TO(i)m(tw)) - (m(t)) (m(^w)) be the 
autocorrelation function for an observable m, R{t,t^) = 
5 {m{t)) / 6h{tv,)\h=f) the linear response oim{t) to a small 
impulse in its conjugate field h at time iw, and x(ij ^w) = 
/j dt' R{t,t') the corresponding response to a field step 
h{t) = /i6(i — iw). In equilibrium, C{t,tvj) = C{t — iw) 
by TTI (similarly for R and x), and the FDT reads 
-^xit - iw) = i?(t, tw) = T al:<^(* - *w), with T the 
thermodynamic temperature (we set fee = 1). A para- 
metric "FDT plot" of X vs. C is thus a straight line of 
slope —l/T. In the ageing case, the violation of FDT can 
be measured by a fiuctuation-dissipation ratio (FDR), 
X(t, t„), defined through 
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x{t,t„) = i?(i,tw) 



T dt^ 



C{t,U) (1) 



In equilibrium, due to TTI, the derivatives -S- in the 
FDT can be replaced by — ^; in the non-equilibrium 
definition (n[) this would make less sense since only the 
t„-derivative of x(i, t^) is directly related to the impulse 
response R{t, t^). Values of X different from unity mark 
a violation of FDT, and can persist even in the limit of 
long times, indicating strongly non-equilibrium behavior 
even though one-time quantities such as entropy and av- 
erage energy may evolve infinitesimally slowly. 

Remarkably, the FDR for several mean field mod- 
els H,^ assumes a special form at long times: Tak- 
ing iw ^ cxD at constant C — C{t,t^), X{t,t^) -^ 
X{C) becomes a (nontrivial) function of the single ar- 
gument C. If the equal-time correlator C(i, t) also ap- 



proaches a constant Cq for i — > cxd, it follows that 
X{t,t^) = J^°^^^^dCX{C)/T. Graphically, this limit- 
ing non-equilibrium FDT relation is obtained by plot- 
ting X vs C for increasingly large times; from the slope 
~X{C)/T of the limit plot, an effective temperature p] 
can be defined as Tefi(C) = T/X{C). In such a scenario, 
where x(^i^w) and C{t,t^) become simple functions of 
each other, either t^, or t can be used as the curve pa- 
rameter for the FDT plot; the latter is the conventional 
prescription [^J^. In general, however, the definition (|^) 
ensures a slope of —X{t, t,„)/T for a parametric x-C plot 
only if tw is used as the parameter, with t being fixed. 

In the most general ageing scenario, a system displays 
dynamics on several characteristic time scales, one of 
which may remain finite as tw — * oo, while the others 
diverge with iw If, due to their different functional de- 
pendence on iwj these time scales become infinitely sep- 
arated as iw — * oOj they form a set of distinct "time 
sectors"; in mean field, Tcft(C) can then be shown to be 
constant within each such sector M. In the short time 
sector {t — tv, = 0(1)), where C{t,t,„) decays from Cq to 
some plateau value, one generically has quasi-equilibrium 
with Teff = T, giving an initial straight line with slope 
-l/T in the FDT plot. The further decay of C (on age- 
ing time scales t—t^ that grow with t^) gives rise to one of 
three characteristic shapes: (i) In models which statically 
show one step replica symmetry breaking (RSB), e.g. the 
spherical p-spin model H , there is only one ageing time 
sector and the FDT plot exhibits a second straight line, 
with Toff > T . (ii) In models of coarsening and do- 
main growth, e.g. the 0{n) model at large n, this second 
straight line is flat, and hence T^s = oo H- (iii) In nrod- 
els with an infinite hierarchy of time sectors (and infinite 
step RSB in the statics, e.g. the SK model) the FDT plot 
is instead a continuous curve Q . 

Tcff has been interpreted as a time scale dependent 
non-equilibrium temperature, and within mean field has 
been shown to display many of the properties associ- 



ated with a thermodynamic temperature g. For exam- 
ple (within a given time sector), it is the reading which 
would be shown by a thermometer tuned to respond on 
that time scale. Furthermore — and of crucial importance 
to its interpretation as a temperature — it is independent 
of the observable m used to construct the FDT plot Q . 

While the above picture is well established in mean 
field, its status in non-mean field models is less obvi- 
ous. To check its validity, one must demonstrate that a) 
a limiting FDT plot exists and that b) it gives effective 
temperatures that are independent of the observable-field 
pair used to calculate C and x- Depending on whether 
there are one or many ageing time sectors one then ex- 
pects FDT plots similar to those for mean field systems. 

Encouragingly, molecular dynamics (MD) and Monte 
Carlo (MC) simulations of binary Lcnnard-Jones mix- 
tures m , as well as MC simulations of frustrated lattice 
gases Q (which loosely model structural glasses, whose 
phenomenology is similar to that of the p-spin model) 
show a limiting plot of type (i). MC simulations of the 
Ising model with conserved and non-conscrvcd order pa- 
rameter in dimension d = 2 and 3 show a plot of type 
(ii) Q . And MC simulations of the Edwards- Anderson 
model in d = 3, 4 ||l^ give a plot of type (iii). The major- 
ity of existing studies, however, do not show that T^e is 
independent of observable (notable exceptions are [0,||) 
since they consider just one observable-field pair. 

In this work, therefore, we analyse a model where FDT 
plots can be calculated analytically for arbitrary observ- 
ables, allowing us to investigate in detail whether the 
simple mean field picture applies. Trap models 11 - 1^] 
are obvious candidates for such a study. Popular as alter- 
natives to the microscopic spin models discussed above, 
they capture ageing within a simplified single particle 
description. The simplest such model |ll[| comprises an 
ensemble of uncoupled particles exploring a spatially un- 
structured landscape of (free) energy traps by thermal 
activation. The tops of the traps are at a common en- 
ergy level and their depths E have a 'prior' distribution 
p{E) {E > 0). A particle in a trap of depth E escapes on 
a time scale t{E) — tq exp{E/T) and hops into another 
trap, the depth of which is drawn at random from p{E). 
The probability, P{E,t), of finding a randomly chosen 
particle in a trap of depth E at time t thus obeys 

{d/dt)P{E,t) = -T-\E)P{E,t)+Y{t)p{E) (2) 

in which the first (second) term on the RHS repre- 
sents hops out of (into) traps of depth E, and Y{t) = 
(r~^(_E))p,^ ^, is the average hopping rate. The solu- 
tion of (H) for initial condition Pa{E) is 

P{E,t) ^ Po{E)e-*^^'^^^ 

+ p{E) /dt'r(t')e-(*-*'^/"(^^ (3) 



For the specific choice of prior distribution p{E) ~ 
exp(— i?/Tg), the model shows a glass transition at a 
temperature Tg. This can be seen as follows. At a tem- 
perature T, the equilibrium Boltzmann state (if it exists) 
is PcqiE) ex t{E)p{E) cx exp(^/r)exp(-^/rg). For 
temperatures T < Tg this is unnormalizable, and cannot 
exist; the lifetime averaged over the prior, (t) , is infi- 
nite. Following a quench to T < Tg, the system never 
reaches a steady state, but instead ages: in the limit 
tw —^ oo, P{E,tv/) is concentrated entirely in traps of 
lifetime r = 0(iw)- The model thus has just one charac- 
teristic time scale, which grows linearly with the age tw 
(In contrast, for T > Tg all relaxation processes occur on 
time scales 0(1).) In what follows, we rescale all energies 
such that Tg = 1, and also set tq = 1. 

To study FDT we extend the model by assigning to 
each trap, in addition to its depth E, a value for an (ar- 
bitrary) observable to; by analogy with spin models we 
refer to m as magnetization. The trap population is then 
characterized by the joint prior distribution a{m\E)p(E), 
where a{m\E) is the distribution of m across traps of 
given fixed energy E. We focus on the non-equilibrium 
dynamics after a quench at t = from T = cxd to T < 1; 
the initial condition is thus Po{E,m) ~ a{m\E)p{E). 
The subsequent evolution is governed by 

{dldt)P{E,m,t)^ 

-T-\E, m)P{E, m, t) + Y{t)p{E)a{m\E) (4) 

where the activation times are modified by a small field 
h conjugate to m as t{E^ m) = t{E) exp (mh/T). Other 
choices of t(E, m) that maintain detailed balance are 
possible [^Q; we adopt this particular one because, 
in the spirit of the unperturbed model, it ensures that 
the jump rate between any two states depends only on 
the initial state, and not the final one. 

The autocorrelation function of m (for h — 0) is 

C{t,tvj) — dm drriyj dE dE^ mm^ P{E^,m^,t^) 

[P{E,ra,t\E^,ra^,t^)~P{E,ra,t)] (5) 

in which P{E,m,t\E^,m^,t^) is the probability that a 
particle with magnetization m^ and energy E^ at time 
tw subsequently has m and E at time t. This obeys 

P{E,m,t\E^,my,,t^) = 

6im - m^)6{E - E^)e-^*-*'^/^^'^^^P{E^, m^, iw) 



dt' 



^(-Bw)e~(*'^*''/"(^''F(-B,i-t')f7(m|£;) (6) 



from which Y{t) has to be determined self-consistently. 



The terms corresponds to the particle not having hopped 
at all since t^ (first term) or having first hopped at t' 
(second term) into another trap; after a hop the parti- 
cle evolves as if "reset" to time zero since it selects its 
new trap from the prior distribution a{m\E)p{E), which 



describes the initial state of the system. We also have 
P{E,m,t) = a{m\E)P{E,t) since at zero field the un- 
perturbed dynamics is recovered. Substituting these re- 
lations into (P) and integrating over m and rrivr, we obtain 
an exact integral expression for C. This is expressible as 
the sum of two components: The first depends only on 
the mean rn{E) of the fixed energy distribution a{m\E) 
and the second only on its variance, A'^{E). 

To find the corresponding response function, we proved 

T^xit. iw) = ^C(i, iw) + ^ (™(i)) (™(iw)) (7) 

in which {m{t)) = ^{E)) p,^ ^^ is the global mean of 
m. Equation (0) generalizes the results of [|lj,0 to the 
case of non-zero means (tti); it is exact for any Markov 
process in which the effect of the field on the transition 
rate between any two states depends only on the initial 
state and not the final one. Substituting our expression 
for C into (M) , and integrating over t^ with the boundary 
condition x(twi^w) = 0, we find a closed expression for 
X- This again comprises two components which depend 
separately on the mean and variance of a(m\E). For 
convenience, we rescale the field h — > Th, thus absorbing 
a factor l/T into the definition of the response function. 
In this way, the slope of the FDT plot of % vs C becomes 
—X = —T/TeS (= — 1 in equilibrium). 

Using our exact expressions, we numerically calcu- 
lated C and X for ^ number of different distributions 
a{m\E), each specified by given functional forms oirn(E) 
and A'^{E). Each form effectively corresponds to a dis- 
tinct physical identity of the observable m. For exam- 
ple a{m\E) = 5{m — E) implies m = E, 'u\ which case 
h — —5T/T. For simplicity, we confined ourselves to dis- 
tributions either of zero mean (but non-zero variance) or 
of zero variance (but non-zero mean). As expected, we 
find that the decay of C (and growth of x) depends not 
only on the time interval t — ty^^ but also explicitly on 
the waiting time t^, consistent with the non-equilibrium 
breakdown of TTI. For our generic observables to, the 
equal-time correlator C{t, t) can also depend strongly on 
t; this is in contrast to the spin models often considered 
in mean field studies, where C{t, t) is automatically nor- 
malized. While from (Q) a plot of x vs C, with t fixed 
and i„ as the parameter, still has slope —X (with t as pa- 
rameter this would no longer be guaranteed a priori)^ the 
variation of C{t,t) will prevent a limit plot from being 
approached as we increase t. We therefore normalize and 
plot X — x/C(t,t) vs C = C/C(t,t); the normalization 
factor C(i, t) is independent of i^ and the same for both 
axes, ensuring that the slope remains —X as desired. 

For the first class of observables (with zero mean 
rn{E)), we considered a variance A'^{E) — exp{En/T) for 
various values of the exponent n, generalizing the results 
of jlj] for n = 0. (For the weaker power law dependence 
A^(_E) = E'^p we can show analytically that the limiting 
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FIG. 1. FDT plots of X vs C for a distribution a{m\E) of 
variance exp{?iE/T) (but zero mean) for n — 0.2, 0.0, —0.2; 
T = 0.3. For each n data are shown for times t — 10^, 10^; 
these are indistinguishable, confirming that the limiting FDT 
plot has been attained. Dashed line: The predicted asymptote 
X = l — Cfori^oo and (7 — > 1. 



FDT plot reduces to that for n — 0.) The equal-time 
correlator is then C(i, t) = JdE P{E, t) exp (nE/T). For 
large i, its scaling behavior can be found using 

/ f^^' expiE/T - E) for t{E) - e^/^ « t 
^ ' '^\e exp{-E) for t{E) - e^/^ » t. 

(which follows from (|), using Y(t) ~ t^-^ ^^)- For 
n < T— I, the integral for C{t, t) converges in a range E ~ 
0(1), giving C{t, t) ^ t'^^^: The equal-time correlator is 
sensitive only to shallow traps (the population of which 
depletes in time as t due to ageing), and we thus 
expect the two-time correlator to decay on time scales 
t—t^ = 0(1), probing only quasi-equilibrium behavior. In 
contrast, for T — 1 < n < T the integrand for C{t, t) has 
most of its weight at energies corresponding to t{E) — 
0{t), yielding C{t,t) ~ t". For such n, the two-time 
correlator will decay on ageing time scales t—tv, — 0(iw), 
and we expect strong violation of equilibrium FDT. (The 
regime n > T is meaningless: It gives C{t,t) = cx)Vt.) 

Our numerical results confirm the above predictions. 
For values of n < T — 1 the normalized FDT plot ap- 
proaches a straight line of equilibrium slope —1 at long 
times t — > oo. In the regime T — 1 < n < T (fig. m), equi- 
librium FDT is strongly violated, as expected; a limiting 
non-equilibrium FDT plot is nevertheless approached at 
long times. This can be proved analytically by showing 
that C and x share the same scaling variable {t—t^)/tvt in 
this limit. The slope of each plot varies continuously with 
C. In contrast to mean field, this is not due to an infinite 
hierarchy of time sectors; the variation is in fact contin- 
uous across the single time sector t — t^ — 0{t^). More 
seriously, different observables give nontrivially different 
plots: at values of C corresponding to a fixed {t—t^)/t^, 
the slopes —X are not independent of n. 

Next we considered a fixed energy distribution a{m\E) 
with zero variance and mean rn{E) = eKp(En/2T). Here 
too, we found that for n < T — 1 the limiting FDT plot 
is of equilibrium form, while for T ~ 1 < n < T (fig. 0) a 
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FIG. 2. FDT plots of x vs C for a distribution with mean 
exp{nE/2T) (but zero variance), for n = 0.2, -0.2; T = 0.3. 
Curves are shown for times t — W^, 10^, but are indistin- 
guishable except for the zoom-inset on the left hand (up- 
per curve: t = lO'^). Dashed line: The predicted asymptote 
X = 1 - C* for i -> oo, C* ^ 1. 



non-equilibrium FDT plot is approached as i — * oo, with 
a shape dependent (now very obviously) on the observ- 
able m. Finally, we considered a power law dependence 
of the mean, rn{E) = E^; p = 1 corresponds to energy- 
temperature FDT. Interestingly, in this case no limiting 
plot exists. This is because the amplitude of the correla- 
tor remains finite as i ^ oo, while the response function, 
for any fixed value of C, diverges as Ini. 

In summary, we have studied FDT in the glass phase 
of the trap model, for a broad class of observables to. 
For observables which indeed probe the ageing regime, 
equilibrium FDT is violated; in most cases, a limiting 
non-equilibrium FDT plot is approached at long times. 
In contrast to mean field, this plot depends strongly on 
the observable. It furthermore has a slope which varies 
continuously across the single time sector i— iw = 0(tw)- 
This shows that in this simple, paradigmatic model of 
glassy dynamics, the mean field concept of an FDT- 
derived effective temperature Toff cannot be applied. One 
may dismiss the trap model as too abstract for this to 
be of general relevance; but by placing traps onto a d- 
dimensional lattice and allowing hops only to neighbour- 
ing traps, a more 'physical' model (of diffusion in the 
presence of disorder) can be obtained. Recent work fl^ ] 
shows the scaling of correlation functions to be unaffected 
by this modification; since (R) would also continue to 
hold, our main results should be qualitatively unchanged. 

Could the idea of a Toff derived from FDT neverthe- 
less be rescued, in this and other non-mean field mod- 
els? Consider first the non-uniqueness of our limiting 
FDT plots. One could argue that in order to probe an 
inherent T^g characterizing the non-equilibrium dynam- 
ics, the statistical properties of the observable must not 
change significantly across the phase space regions vis- 
ited during ageing. (In coarsening nrodels, similar argu- 
ments have been used to exclude observables correlated 
with the order parameter H.) Since in the trap model 
the typical trap depth E increases without bound for 



t ^ oo, a 'neutral' observable would presumably require 
that rrT{E), A'^{E) — » const, as i? — >• cxd; with this restric- 
tion we indeed find a unique limiting FDT plot. Its slope 
~X will still vary continuously with t — t^f,, however, ex- 
cluding the link to a thermodynamically meaningful Toff. 
Similarly 'rounded' FDT plots have recently been found 
in coarsening models at criticality Pq] ; the limiting value 
—Xao of the slope for C ^ was there shown to be a uni- 
versal amplitude ratio. It is possible that at least this X^o 
could define a sensible Toff, and in fact all our limiting 
FDT plots above share a common value X^o = 0. In con- 
clusion, there is at least one class of 'critical-like' models 
(which includes the trap model |l7[ ) for which a non- 
equilibrium Toff can be defined at best from the asymp- 
totic FDT Xoo, and possibly only for sufficiently neutral 
observables. It remains an open challenge to character- 
ize this class more fully, and to delineate it from those 
models which appear to exhibit mean field-like behavior 
in spite of their non-mean field nature (7|-|9|] . 
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